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Abstract 

We use computational algorithms recently developed by us to study 
completely four index divergence free quadratic in Riemann tensor polyno- 
mials in GR. Some results are new and some other reproduce and /or cor- 
rect known ones. The algorithms are part of a Mathematica package called 
Tools of Tensor Calculus (TTC)[web address: http: / /baldufa.upc.es/ttc 



1 Introduction 

Recently the interest of tensors containing the interaction helds, and in particu- 
lar the Riemann tensor, has been increased. The so-called superenergy tensors 
are good examples, having some particular properties [EJ. This interest is not 
new. The tensors of Bel- Robinson Bel || and Sachs B among others are some 
examples which have been studied in the past and having zero divergence in 
certain conditions. CD Collinson has made a much deeper study giving all 
four index quadratic in Riemann and divergence free tensors in four dimension 
for any space-time. Nevertheless this kind of calculations are very hard to do 
by hand and some errors could arise. Having this in mind, recently, we have de- 
veloped some computational algorithms able to handle index tensor polynomial 
expressions. The paradigm of these polynomials are those built with Riemann 
tensor as in the example 

_J_o mn p kl ■ o p fcmp / 
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The implementation of the algorithms are part of a package called TTC || 
( Tools of Tensor Calculus) . These algorithms allows us to perform a complete 
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simplification of this kind of expressions. The simplification algorithms do not 
find the most beautiful expression but, besides zero result, a canonical one. 

Other authors who have worked on algorithmic simplification of index tensor 
polynomials are Ilyin and KryuchovQ, Portugal^, Fulling et al.|J. Parker and 



Christiensenjll do not fall on this list because although they have a big pack- 
age on tensor calculus, they don't work algorithmically in this topic but with 
libraries of rules. In a recent article[[lOj we have solved the part of the simplifi- 
cation problem related to dummy indices and monoterm symmetry properties. 
Cyclic and in general multiterm symmetries properties were not considered. 
Also dimensional properties were left aside for further works, due to its com- 
plexity. Actually [l| we can handle multiterm and dimensional properties 
with the advantage that the algorithm includes the possibility to store a com- 
pacted version of the generated rules in order to increase functionality. This 
means that, although a first computation over a polynomial can be really slow, 
others on the same or on other sessions are faster than the first one. This allows 
us to play easily with some results by permuting indices or making some other 
operations. 

There are some trust problems when using programs which do automatic 
computations, specially when no exact coindicences arise when compared with 
older papers. All we can say is that our algorithm has been tested, besides 
in the present paper, on a large set of known Riemann properties ( among 
others all those included in and calculations related to the Gauss-Bonnet 
invariant . On the other hand the fact that the implementation of the algo- 
rithm is part of a package, TTC,{ from now on we will use the name TTC to 
call both, the simplification algorithm and/or the package) it allows us to test 
the results on explicit space-times. We have tested the calculations on Vaydia 
metric confirming our results. 

The aim of this paper is: 

• To show the TTC capabilities on a real calculus. This can be of inter- 
est to researchers working on algebraic computation and/or working on 
superenergy tensors or on any other field where these techniques can be 
applied. 

• To review some known results and correct some others. Specially we have 
corrected a minor part of the results of CD. Collinson and some sign on 
the original Sachs paper 

• To prove some theorems about existence and uniquness of superenergy 
tensors. 

By superenergy tensors we understand, in this paper, four index quadratic 
in Riemann tensor having some chosen properties. Following the work of CD. 
Collinson we will find divergence free tensors and from these we will analyze 
the index symmetry properties. We leave aside other properties such as the 
dominant super-energy property (DSEP) and others explained in reference (T) 

The plan is as follows. In section 2 we briefly explain what TTC does when 
simplifying index tensor polynomials in general. Also we describe some TTC 
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functions usually used in operations over this kind of objects . In section 3 
we define exactly what kind of tensors we will find and which are the steps 
to reach our main result, that is to find all four index quadratic in Riemann 
and divergence free tensors in any dimension and in any space-time. Next we 
specialize these results for four dimensions and for vacuum space-times. We 
analyze the tensors found with respect to index symmetry properties, specially 
we give existence and uniquess theorems for totally symmetric tensors. 

Finally Appendix A is devoted to show a sample TTC session where the 
C.D.Collinson calculations are reproduced using the essential TTC functions. 



2 Computational simplification of index tensor 
polynomials 

As stated above, by simplification we mean to find a canonical, therefore unique, 
index tensor polynomial for any given monomial. In this process the following 
must be taken into account: 



Dummy indices: the output expression must be unique with respect to 
the freedom of renaming dummy indices. 

Index symmetry properties of each singular tensor appearing in the mono- 
mial: since in this paper we work with Riemann tensor R the following 
properties, and those which arise when diferentiating and contracting, 
must be taken into account 

Rijkl = —Rjikl = —Rijlk = Rklij 
Rijkl + Riljk + Riklj = 
Rijklym ~t~ Rijmk:l ~t~ Rijlm\k 

= a Rijki- n....-m +(Riemann terms) .... (1) 

k k-2 

where the last one refers to Ricci commuting relations. At this point it is 
important to set up the convention for Ricci commuting relations and the 
definition of Ricci tensor used by TTC which follows [[15) . 

If v is a vector field then 

v i;j\k — v i;k:j = V m -fT™ j k 

p _ pm 

ij imj 

It is important to note the change of sign of the Ricci tensor with respect 
to that used in Sachs @, Collinson§ and Robinson [fl7| original papers 

Properties due to the finite dimension of space. This kind of properties 
arise due to the posssibility of antisymetrization over n + 1 indices appear- 
ing in the monomial, n being the dimension of space. In the next example 
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ii = 2 

... + A i:j B k + ... => A [Lj B k] =0 

In the following subsection we briefly explain how algorithms work. For details 
we refer to |l§ [l3| 

2.1 Dummy index problem 

Actually TTC uses a codification of each monomial, which is an improved ver- 
sion of |iC| |, following the criteria: 

1) Free indices are renamed: IndexF [1] , IndexF [2] , ... all of them in up- 
per position ( contravariant) 

2) Dummy indices are renamed IndexD[l] , IndexD[2] , . . . all of them 
in upper position ( contravariant) 

3) We find the permutation of free and dummy indices in order to find the 
most lexicographically ordered monomial expression. Note that at this stage 
the codification is independent of the order of free indices. We call this global- 
codification 

Finally we take the invers of permutation performed on free indices over the 
globalcodification of the monomial. The final result is called codification of the 
original monomial. 

All rules found through a calculation are stored in codification version but 
all monomials with the same globalcodification version can use a permutation 
version of the stored rules, so the work of finding rules is essentially independent 
of the order of free indices. 

2.2 Index tensor properties 

First the user must declare each tensor entering the monomial, together with 
their symetries. Some properties are automatically applied as Ricci commuting 
relations. 

Key functions: InputTensor, InputSymmetries 

TTC acts on each monomial applying in all possible ways the individual 
properties of each tensor over the starting monomial and over each new resulting 
monomial. When this loop ends ( no new monomials are created) TTC solves the 
resulting system of equations and applies the rules over the original monomial. 
Subsequent monomials can use the same or a permuted version of these rules if 
the globalcodification version is the same as one of the monomials in the rules. 

Key functions: Index [gname , Simplif yAlllndex [2] ] 

2.3 Dimensional properties 

As stated in the introduction, dimensional propeties arise due to the possibility 
of antisymmetrizing n + 1 indices of the starting monomial and equating the 
result to zero, n being the space dimension. Nevertheless we must be sure to 
take all possible sets of n + 1 indices in the monomial, even if they are dummy or 
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hidden ( as in the curvature R). The algorithm takes the monomial and writes 
it in a metric expanded version. The possible sets of n + 1 indices are in the 
metric part of the monomial. Let us see some examples in four dimensions. 

1) In the first example there are no dimensional properties: 

-> R Q mp Q nq 

2) Only one set of 5 indices can be effectively antisymmetrized: 

p pab p p n mn n ar r,bs oq pr 

titi > ti mrns K opqr g g g g H g F 

3) The last example is the case handled in the present paper. Several sets of 
5 indices can be efectively antisymmetrized. Note that for n > 6 there are no 
dimensional properties. 

titi — > ti mnop tt rstv g g y g g g g 

TTC finds all these properties and simplify them using nondimensional ones. 
Then it solves the resulting system of equations and applies the resulting rules 
when necessary. Although the procedure described is correct actually the algo- 
rithm has been improved in order to optimize the number of sets of indices to 
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Key functions: Index [gname , Simplif yAlllndex [3] ] 



3 Collinson like calculus 

CD Collinson found by hand 10 independent divergence free tensors having four 
indices and quadratic in Riemann. Here we use TTC to make the computerized 
version of his calculations, but the plan is close to the Collinson paper. 
First we must generate the starting basis 

rpabed Aabcdefghklmn p p i ryabcdklmnpq p 

^ Si -R-efgh ^klran ~T~ ' J FLklmn-pq 

where A abcde fs hklmn and B abcdklmn Pi are polynomic tensors built with the met- 
ric g ab using the free indices plus the contracted ones in all possibles ways, 
symbollically 

j^abcdefghklmn = ^ a ^ a (^gib gCd f ^gh gkl gW,n\ 
a 

fiabcdklrnnpq = f3 a a (g ab g cd g U J m " g pq } 
a 

a being the permutations over metric indices and a a and /3 CT constant parame- 
ters. 

TTC performs this calculation simplifying the result using nondimensional 
properties and dimensional ones. In the second case we obtain, as in Collinson's 
paper, 63 independent monomials. This polynomial is what we call the starting 
basis T, which has 63 independents coeficients (including the global one). In 
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doing this TTC finds 3 independent dimensional properties which includes those 
found by H.A.Buchdahl|]l|. 

Second we take the divergence of T 

rpabed 

;a 

TTC calculates this and simplifies it by using nondimensional and dimensional 
properties. After this we must solve the coefficients of the equation 

T abcd . a =0 (2) 

so we find the family of polynomials T fulfilling (j^) . This computation has been 
done in a generic space-time and in vacuum. The results have been analyzed 
under a few relevant index symmetry properties. From the computations we 
have: 



3.1 Without-dimension properties 

Theorem 1 In a generic space-time and not using dimensional properties: 

1) There exist 14 independent quadratic in Riemann four index divergence 
free tensors. 

2) There are no totally symmetric tensors fulfilling 1 ) 

3) The complete family of tensors T abcd fulfilling 1 ) totally symmetric in 
(bed) is 

Tabcd ~ rpabed i ~ rpabed. 
— a s l s +a R l R , 

rpabed Qa(bcd) . 

Qabcd \ gac R d R ib _j_ 2 R ab:cd 4 R bd;ac_^_ 

\g ac R bd ' 1 i - 2 g ac R b + § R b t R acld + 

2 R a ^3 R c d — g ac R. d _ ffijbk 

Tiabcd vab ~cd. 
r — A 9 ) 

X ab — K U ab + L V ab - \W ab 

rjab =G ab ' s s ~2G sb ' a s + 2G a p Rv b -\G pq R pq g ab 

yab = R ;ab _ R :s ^ g ab _ R gab 

w ab = n a Pi r u b _ I n ab n m P ar r 

r<ab _ r>ab a [° q b \ 

° cd~ K cd 4 3 [c b d] 
gab = R ab _ l g ab R 

where as, an, K, L are four independent constant parameters. 

Ts is the tensor found by Sachs Q with an opposite sign on the term § R b l R ac 
with respect to his original definition. Tr is the family defined by Robinson JlT[| . 

Theorem 2 In a vacuum Rij = space-time and not using dimensional prop- 
erties: 



G 



1) There exist 6 independent quadratic in Riemann four index divergence 
free tensors. 

2) There is only one tensor T abcd fulfilling 1 ) and totally symmetric: 

rpabed j^(abcd) _|_ ^g pa(bcd) 4 Q{abc)d 4 Q(a\d\bc) . 

M abcd = g ab g cd R mnop R^nop. 



abed p a b TDmcnd. 

abed „ab p c ryranod 



Qabcd „ab p c pn 

3.2 Using dimensional properties 

Theorem 3 (Collinson corrected) In a generic four dimensional space-time: 
1) There exist 9 independent quadratic in Riemann four index divergence 
free tensors. 

Using the same notation of Collinson ^/ the Collinson family is 

Tabcd ~ rpabed i ~ rpabed i ~ rpabed i ~ rpabed i ~ rpabed , ~ rpabed i 
— a\l 1 +a2± 2 +03^3 +a4l 4 +ac,l 5 +aQi 6 + 

„ rpabed i ~ rpabed i ~ rpabed 

agJ 8 +a9ig +aioi 10 

where etj are arbitrary numerical coefficients. 

The tensors T\ to Tg ( without T-? ) can be writen in terms of T 10 



1 n cdrpa i b 



and Tio is defined through Tf 6cd = 6 Qf cd + Qf (cd) + Qf (cd) being 

Qabcd _ j-yc j-ybdai j^b j^daci _ j-yaibc \d _|_ j^acdi \b _|_ g ab j-ydicj 
Qabcd = 4 ^aibd c _ g ^ad^b ;c. _ g g ad R icbj^ + g fiai ,; cd + g R abcijd 
Qabcd g q ac Tl ^ J^jkbi g^£ cc ^ J^biaj _|_ g^£ d j^bica _ g^j> 6d j^cia j j 



rpabed 


rpacbd 
1 2 — 


rpadeb 
1 3 — 


r-cdrpai b 

9 1 io i 


rpabed 
1 6 — 


rpacbd 
1 h — 


rpadeb 
1 A — 


n cdrpa b i 
9 i 10 i 


^ rpabed _ 


rpacbd 
~ A \Q 


rpabed 
~ 1 10 




rpabed 
1 9 — 


rpadeb 
-^10 


o rpabed 





2i? d , ( i?^ c - 2i?/ b j i? cia ^' - Sg^R^R" + 3g ab 'R^l R jkc 



2) There are no totally symmetric tensors fulfilling 1 ) 

We note that the Collinson tensor T7 does not appear, which is identically 
zero, since Tiq is in fact symmetric with respect to the last two indices (this is 
not mentioned in Collinson's paper). 

Theorem 4 In a vacuum Rij = four dimensional space-time: 

1) There exists 1 quadratic in Riemann four index divergence free tensor. 

2) The tensor found in 1) is totally symmetric and is the Bel-Robinson 
tensor ||/. 
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4 Concluding remarks 



Trough four theorems we have reviewed and improved some results on classical 
gravitational superenergy tensors. The results have been reached using compu- 
tational algorithms which do automatic calculations. As we have commented in 
the introduction some trust problems could arise. 

Nevertheless we think that trust like problems are not excusively related 
to computational techniques. How can we be sure that all calculation of CD. 
Collinson, or anyone, are correct ?. Collinson is a good pioneer on large calcula- 
tions on the topic but to validate his results we hope some other should obtain 
the same results using independent calculations. We also hope that any other 
program will be able to do our calculations with an independent code and/or 
algorithms... 

Assumming that these computational techniques are validated and limiting 
our scope to superenergy tensors we think that they are a very good comple- 
ment to research based on defining such tensors having some choosen properties. 
We can analyze the uniqueness of these definitions with respect to the desired 
properties. But we can also analize more complex problems, specially those 
related to the existence of conserved quantities. These studies can be done on 
a generic space- time and using Einstein field equations. We hope that in the 
near future we will be able to handle the case of cosmological vacuum and scalar 
field interaction. 
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Appendix: The TTC master session 

This appendix is devoted to show a sample of TTC session so that somebody 
else can reproduce the calculations using TTC. In this session we find the result 
of Theorem^). The inputs are indicated by In[] :=. The outputs have no 
special sign. The comments are enclosed by (* comments *). 

(* first we load the program ttc.m. we load also eines.m 

where some private tools are defined.*) 

In[]: 

<<ttc .m 

<<eines .m 

I Tools of Tensor Calculus 4.1.0 I 
I A.Balfagon, P.Castellvi and X.Jaen I 

I \protect\vrule widthOpt\protect\href {http : //balduf a. upc . es/ttc}{http : //balduf a. upc . es/ttc} I 

I e-mail: ttcObaldufa.upc.es I 
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version: September , XX, 1999 I 



I Session started on I 

I September, XX, 1999 I 

I at 16 h 8 min 52 s I 



CPU Time=0. s Memory in use=2.434 Mb 



Eines Personals 1.0 I 
X.Jaen & A.Balfagon I 
version: juny, 23, 1999 I 



(* Here we declare the name of the coordinates and dimension to be used when needed. We also introduce a 
metric and a Riemann tensor. 

Finally we introduce the indices we want to see in outputs *) 
In[] : = ( 

InputCoordinates [cx4,4] ; 
ScalarBasis=cx4 ; ScalarBasisQ=False ; 
InputSMetric[gn,cx4,"g",g] ; 
InputSRiemann[gn,cx4, "R" ,Rie,Ric,R] ; 
Input Index [{i , j , k, 1 ,m,n, o,p,q}] ) ; 

(*This is usefull if we want to store in "filename" all rules 
generated by TTC. If you use it remember end save using 
EndSimplif yAlllndexSave*) 

(*IniSimplif yAlllndexSave [filename] *) 

(* generator is where the kind of tensors we study are defined*) 
In[] :=generator= 

Rie [a,b, c,d] Rie [-a,-b,e,f ] +Rie [a, c ,d, e , . ; -a, . ;f ] //Index [gn] 

ijkl kl mnij mijk .;1 

+ R R + R 

m n . ;m 

(* IndexTensorBasis takes the generator and uses some easy tricks 
to build the starting basis, simplifying them using 

nondimensional and dimensional properties. It uses A [number] 

as parameters. The result is stored as base*) 

In [] : =generator//Index [gn, IndexTensorBasis [3, A, base , CoRiemannRules [gn, cx4] ] ] ; 
(* see *) 

(*EndSimplif yAlllndexSave*) 

(* we can see the starting basis) 

In[] :=base 

ijkl .;k .;1 i j .;j .;1 i k 
+ A[l] R g + A[2] R g + 

ikjl mnop ijkl mnop 

> A[62] g g R R + A[63] g g R R 

mnop mnop 



(* Here we define A[l].. A[63] as constants parameters*) 
In [] : =SetConstantCoef f [A , cx4 , 63] ; 

(* Here we compute the divergence of base and simplify up to 
dimensional properties assigning the result to dbase *) 

In[] : =dbase=base [ . ; -i] / /Index [gn, Simplif yAlllndex [3] ] ; 

(* Here we find the parameters which solves dbase=0 assigning 
the result to rules*) 
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In[] :=rules=SolveCoeff [dbase.A] [[1]] ; 



(*Here we apply rules to the starting basis base assigning 
the result to baseXA*) 

In[] :=baseXA=base/. rules; 

(* baseXA fulfill the condition imposed...*) 

In [] : =baseXA [ . ; -i] //Index [gn , Simplif yAlllndex [3] ] 



(* Here we find tlxa . . .t9xa which are the nine independent tensors 
we have found corresponding to the nine independent coefficients A [number] 
appearing in baseXA*) 

In[] : =tlxa=baseXA/ . 
{A [6] ->1 , A [9] ->0 , A [12] ->0 , A [58] ->0 , A [59] ->0 , 
A [60] ->0 , A [61] ->0 , A [62] ->0 , A [63] ->0} ; 

In[] :=t2xa=baseXA/. 
{A [6] ->0 , A [9] ->1 , A [12] ->0 , A [58] ->0 , A [59] ->0 , 
A [60] ->0 , A [61] ->0 , A [62] ->0 , A [63] ->0} ; 

In[] : =t3xa=baseXA/ . 

{A [6] ->0 , A [9] ->0 , A [12] ->1 , A [58] ->0 , A [59] ->0 , 
A [60] ->0 , A [61] ->0 , A [62] ->0 , A [63] ->0} ; 

(* also for tx4,tx5,tx6,tx7,tx8*) 

In[] : =t9xa=baseXA/ . 

{A [6] ->0 , A [9] ->0 , A [12] ->0 , A [58] ->0 , A [59] ->0 , 
A [60] ->0 , A [61] ->0 , A [62] ->0 , A [63] ->1} ; 

(* we can see some of them*) 
In [] : =t4xa 

ijkl 2 i 1 j k 2 i k j 1 jlik 

+ Rg g - Rg g +2 Rg R 

jkil iljk ikjl jl imk 

2 Rg R -2 Rg R +2 Rg R -2 g R R 

m 

il jmk jk iml ik jml 

4 gRR+2 gRR-4 g R R - 

m m m 

iljk mn ikjl mn ik.;j.;l 

2 ggRR+2 ggRR-2 R + 

m n m n 

i 1 . ;j . ;k j k . ;i . ;1 iljk . ;m 

2R +2R -2gR 

. ;m 

j 1 .;i .;k i k j 1 . ;m 1 m i j k 

2R +2gR -2RR + 

. ;m m 

kmijl jmikl ilmjnk 

2RR +2RR +2R R 

m m m n 

ikmjnl ilmnjk ikmnjl 

2R R ' - R R +R R + 

m n m n m n 
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1 J 



m n k 1 



m n o k 



(* now we supose that we have find baseC which is 

the baseXA expressed ussing the basis found by Collinson. 

We use B[l] . . (B[7]) . . .B[10] to name the parameters ( without B[7]) 

accordingly to Collinson notation. We take baseC and use 

IndexUpdate to define T in such a manner that we can use T 

as a tensor in their own right. This will permit us 

to study symmetry index*) 

In [] : =IndexUpdate [T, " :=" , baseC] 

(*Here we study the symmetry between the first two indices, 
simplifying them up to dimensional properties and finding 
the condition for the parameters B[n]*) 

In[] :=T[a,b,c,d]-T[b,a,c,d]//Index[gn,SimplifyAllIndex[3]] ; 
In[] :=rules=SolveCoeff re.B] [[1]] 

{B[4] -> 0, B[5] -> 0, B[2] -> 0, B [3] -> 0, B [8] -> 0, B [9] -> 0, B[10] -> 0} 

(* Only B[l] and B [6] are independent. 
It is in this way how we have analized 
symmetry properties in this paper*) 
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